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Abstract. We have investigated metallic break junctions of the heavy-fermion 
compound UEM2AI3 at low temperatures between 0.1 K and 9K and in magnetic fields 
up to 8T. Both the current-voltage /(V) characteristics and the dV/dI(V) spectra 
clearly showed the superconducting (T c ~ 1.8 K) as well as the antiferromagnetic 
(Tn —14 K) transition at low temperatures when the bias voltage is raised. The 
junctions with lateral size of order 200 nm had huge critical current densities around 
5 x 10 10 A/m 2 at the antiferromagnetic transition and hysteretic I(V) characteristics. 
Degrading the quality of the contacts by in situ increasing the local residual resistivity 
reduced the hysteresis. We show that those hysteretic I(V) curves can be reproduced 
theoretically by assuming the constriction to be in the thermal regime. It turns out that 
these point contacts represent non-linear devices with A-shaped I(V) characteristics 
that have a negative differential resistance like an Esaki tunnel diode. 
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1. Introduction 

Point-contact (PC) spectroscopy is widely used to study the interaction of conduction 
electrons with elementary excitations or quasiparticles in conducting solids [TJ |2J- 
Energy- resolved PC spectroscopy is possible when the inelastic relaxation length of 



electrons in the contact region A = min{/ in , J / e i/i n /3} (here l e \ and /i n are the elastic 
and the inelastic mean free path of the electrons) is larger than the size or diameter 
d of the contact. In the opposite case of A <C d and a short phonon scattering length 
Z p h < d, the excess electron energy dissipates in the constriction. This Joule heating 
increases the temperature inside the contact when a bias voltage is applied |S1 IH E] • 

Therefore the interpretation of the PC data requires to find out the regime of charge 
transport. l e \ does not depend on energy and can be determined rather accurately for 
the PC region. / in depends on energy, and no method exists to calculate it reliably. 
To identify the transport regime becomes especially important for PCs with complex 
systems like heavy-fermion, high-T c , or Kondo-lattice compounds that typically have 
large electrical resistivities because of their strong electron correlations. 

A. Wexler [Hj derived 

for the PC resistance R as function of temperature T and contact size d. The parameter 
(3 ~ 1 varies slowly as function of l e \/d, and (3 = 1 for large contacts d 3> l e i- Wexler's 
formula interpolates between the ballistic Sharvin (I = l e \ 3> d) and the diffusive Maxwell 
(/ <C d) resistance. The latter describes transport as in the bulk material. 

Since Sharvin's resistance does not depend on temperature, differentiating Eq. (Q) 
with respect to temperature yields 

< 2 » 

for the size of the contact. This is considerably more reliable for deriving d than Eq. Q 
itself. The main reason is that the residual resistivity in the PC region can strongly 
differ from the bulk po, for example due to the stress exerted while forming the contact. 
Eq. (0) was experimentally verified for PCs with simple metals by Akimenko et al [7 . 

The same method can be applied to heavy-fermion compounds. They show at low 
temperatures power-law dependencies of their electrical resistivities p{T) = p + AT 11 
(n = 1,2,3 for the various compounds investigated) which was also revealed in the PC 
resistances |8j. Heavy-fermion compounds typically have large A coefficients because of 
their strong electron correlations, which makes it straightforward to measure dR/dT just 
above T c . For those high-resistivity superconducting (SC) metals the local normal-state 
residual resistivity in the PC region 

p = d5R (3) 

can be extracted from the drop SR of the contact resistance due to SC. Such a 
relationship has been found for a number of heavy-fermion SCs over a wide range of 
contact sizes P ITIH ITT|. 
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Figure 1. Electrical resistivity p(T) of polycrystalline UPCI2AI3 (solid curve an d 
of two single crystals along the basal ab-plane (dotted curve E], dashed curve 
|15p. The inset shows p{T) at low temperatures. Arrows mark the SC and the AFM 
transition, respectively. 



Here we present experiments on PCs between two pieces of the heavy-fermion 
compound UPd 2 Al 3 ^2]; using the technique of mechanically-controllable break 
junctions. Compared to the conventional spear-anvil type technique to form point 
contacts, break junctions have much better mechanical stability. But more importantly 
breaking the sample at low temperatures in the ultra-high vacuum region of the 
refrigerator avoids oxidation of the freshly broken surfaces of the contact interface. 
UPd 2 Al 3 becomes antiferromagnetic (AFM) at T N ~14K. It is SC below T c ~ 1.8 K. 
We have observed huge non-linearities of the PC resistances and even hysteretic I(V) 
characteristics. We derived the contact size and the residual resistivity in the PC 
region according to Eqs. (j2j) and Q, respectively. We found that the very short elastic 
mean-free path in the constriction l e \ <C d points to at least the diffusive regime of 
electron transport through the PC. Considering also the small inelastic mean-free path in 
UPd2Al3, reflected by the steep p{T) rise with temperature around the AFM transition 
in Fig. [TJ we applied the thermal model developed in Refs. jU |3] for the case / e i, l in <C d 
to take into account the locally increased temperature in the PC when a bias voltage 
is applied. Using the experimental p{T) in Fig. [TJ this model described quite well the 
observed I{V) characteristics and their modification with temperature, reproducing also 
the hysteretic features. 

2. Experiment 

We have investigated three UPd2Al3 single crystals. Two of them had one long side in 
the c-direction of the hexagonal crystal lattice, one sample had it in the perpendicular 
ab-direction. A 0.5 — 0.7 mm deep notch was cut into the middle of the ~ 1 x 1 x 5 
mm 3 large UPd2Al3 bars using a diamond saw. This defined the break position. Each 
sample was glued electrically isolated onto a flexible metal bending beam. Twisted pairs 
of voltage and current leads were attached with silver epoxy to both sides of the sample, 
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Figure 2. Resistance R of UPCI2AI3 break junctions along the ab-direction versus 
temperature T before breaking (bottom curve) and with increasing PC resistance. The 
curves, except the upper one, are scaled along the i?-axis to fit into the same window. 
The resistance R n at T > T c just above the SC transition is indicated for each curve. 
The contact of the upper curve is about 200 nm wide (see text). Dashed horizontal 
lines indicate missing data between 2.2 and 3K in some of the R{T) curves. This was 
due to an instability of the mixing chamber of the refrigerator while slowly warming 
up. The inset shows R(T) vs T 2 of the unbroken sample. The straight dotted line 
describes the contribution of the A coefficient. The current excitation was 1 = 1 mA 
for the contacts with R n < 5nrf2, 0.5 mA for the 14mi7, and 2.5 for the 0.66 
contact. It was chosen small enough to not degrade R{T). 

which was then mounted onto the mixing chamber inside the vacuum can of the dilution 
refrigerator. The temperature could be varied between 0.1 and 9K. With a micrometer 
screw the bending beam is bent at low temperatures, breaking the sample at the notch. 
The resistance of the break junction, that is its lateral contact size, could be adjusted 
mechanically with the micrometer screw. For further details of the experimental setup 
see Refs. HH1HI]. 

At room temperature the resistance of the samples with the notch was about 
5 mf2, corresponding to the approximate geometrical cross-section and a contact size of 
0.2 mm. Note that the notch only defines the macroscopic position of the junction, the 
microscopic contact is less well defined. After removing the sample from the refrigerator, 
the surface of the junction was not mirror-like or smooth as expected for a single crystal. 
The fracture was usually tilted with respect to the direction of the notch, and thus the 
crystal axis. Therefore the current flow through the contact might deviate slightly from 
the direction defined by the long side of the sample. Magnetic fields up to 8T could 
be applied perpendicular to the bending mean, that is perpendicular to the long side of 
the samples and to the ideal direction of current flow. 

The I(V) characteristic and the differential resistance dV/dI(V) were recorded by 
injecting a DC current I superposed by a small alternating current 51 and measuring 
the differential voltage drop V. Its alternating part 5V was detected using the standard 
lock-in technique. 
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3. Results 

All three UPCI2AI3 single crystals showed qualitatively the same results. Therefore we 
concentrate here on one of them, that with the long side in ab-direction. Figure El shows 
the temperature dependence of the resistance R(T) of the break junctions below 9 K 
before breaking and of several contacts after successively reducing the contact size by 
increasing the bending force. The superconducting transition at 1.8 K as well as the 
~ T 2 increase above T c like in the bulk samples is clearly seen. Occasionally R(T) 
changes in small steps. The reason for this is that UPd2Al 3 single crystals are quite 
brittle. They also have a large thermal expansion with respect to the bending beam 
above ~ 1 K. When the temperature changes, the stress in the contact region changes. 
This stress is sometimes partly released, slightly varying the contact size or the local 
residual resistivity and, thus, R(T). 

With increasing PC resistance the SC transition broadens. We believe that this is 
mainly due to the stress in the PC area when the sample is broken and the contact being 
formed. Additional broadening is caused by the extremely small critical supercurrent 
which suppresses the Sharvin resistance at low temperatures and small excitation 
voltages, see also the discussion below. On increasing the temperature the critical 
current decreases, so that Sharvin's resistance is again added to the total resistance. 
However, its contribution to R(T) is small since for the investigated contacts Sharvin's 
resistance is much smaller than Maxwell's resistance. Changing the force on the bending 
beam changes the contact size and the stress there, too. This allows us, although in an 
uncontrolled manner, to vary in situ the local resistivity at the PC. 

FigureElshows for the UPd2Al 3 break junction with R n = 0.66 Q as example how the 
I(V) curves typically change with temperature. At low temperatures /(V) is strongly 
hysteretic. At higher temperatures the hysteresis smears out and transforms into an 
inflection point that corresponds to the pronounced dV/ dl maxima above about 5 K in 
Fig. Ufa). Large magnetic fields up to B = 8 T only slightly modified the /(V) curves at 
0.1K by reducing the size of the hysteretic loop. A 4T field as well as the temperature 
above T c completely suppressed the superconducting features, a zero-bias minimum of 
the differential resistance accompanied by a series of spikes, see Fig. Efb). 

4. Discussion 

We start the analysis by deriving the size d of the contacts. Above T c both the specific 
resistivity and the contact resistance vary with the same A T 2 power law. According to 
Eq. (J2J) the contact size d = A bu i k /A PC . Literature values for the A buU , coefficient range 
from 0.15 to 0.25 fiQcmK~ 2 , for example in Refs. ^21 El E3- In part this variation 
could be due to micro-cracks in the bulk samples which spoils the geometrical factor. 
Therefore we choose the average A bu ik = 0.20 /ificmK -2 , which coincides with that in 
|14j . The absolute error in d can then amount up to about ±33%, but the relative 
accuracy needed to compare the different contacts is much better. In this way the 
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Figure 3. I(V) characteristics of the UPC12AI3 break junction with R n = 0.66 at 
the indicated temperatures. Solid (dashed) lines correspond to sweeps with increasing 
(decreasing) current. The hysteretic loops become smaller when the temperature rises, 
vanishing above ~ 5 K. The upper inset shows /(V) below (0.1 K, solid line) and above 
(2K, dashed line) the SC transition in extended scale. The lower inset shows part of 
the I(V) curves at T=0.1K and the indicated magnetic fields. 



contact in Fig. |3] has d ~ 200 nm. 

We can now directly read off the critical current density from the I(V) data in 
Fig. El For the AFM transition the current density reaches up to 5 x 10 10 Am -2 . At the 
SC transition, marked by the dV/dl maximum in Fig. IHb), the critical current density 
approaches 1.5 x 10 10 Am -2 . Both values are lower bounds for the kinetic critical current 
densities because they include local heating of the PC discussed below. 

According to Eq. (j3J) the 5R = 0.66 Q resistance drop due to SC results then in a 
normal-state residual resistivity po ~ 13 /iOcm. This is about three times larger than 
the bulk po, estimated from R(T) of the unbroken junction in Fig. |21 

This 5R includes a possible contribution from the Josephson effect: The differential 
resistance vanishes completely within a very narrow (~ 10 /xV) voltage range around 
zero bias, barely seen in Fig. 0] The much broader (~ 0.3 mV) zero-bias minimum 
has a plateau of around 0.10 Q, fitting well the ballistic Sharvin resistance calculated 
using the known contact diameter. This agreement supports our interpretation that we 
are dealing not with multiply connected contacts but with single contacts. Taking into 
account Sharvin's resistance would slightly reduce the calculated local residual resistivity 
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Figure 4. (a) Differential resistance dV/dl of the UPcb AI3 break junction from Fig.[3| 
at 0.1, 3, and 7K at high biases. Below ~ 5K some of the curves are discontinuous 
around ±5 mV, indicated by the dashed lines, (b) dV/dl of the same break junction 
at low biases. The SC anomaly has disappeared at 3 K or at £?=4T, that is well above 
either T c -1.8 K or B c -3.5 T of UPd 2 Al 3 . 

from 13 fiQcm to 11 pQcm. 

The elastic electron mean free path at low temperatures can be estimated using 
the typical metallic pi ~ 2.5 x 10~ 15 Qm 2 (here I is the elastic mean free path 
and p and I values taken from |12j ) as l e \ ~ 20 nm. This leads to the inequality 
l e i ~ 20nm^C 200 nm~ d for two of the important length scales of the constriction, 
implying that these PCs are at least in the diffusive regime. However, heavy-fermion 
compounds typically have a large residual resistivity and/or already at low temperatures 
a strongly increasing electrical resistivity, they are very likely in the thermal regime |16|. 

The ballistic Sharvin resistance is then negligible, and the PC resistance can be 
described by Maxwell's resistance 

R{T) ~ p(T)/d . (4) 

In contrast to a ballistic contact, where energy dissipates far away form the contact 
region, now all energy is released in the constriction. This increases its temperature 
with bias voltage. Assuming the Wiedemann-Franz law to be valid, the temperature in 
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Figure 5. (a) Smoothed p(T) (symbols) used for modelling the PC. The solid 
line shows Geibel's p(T) from Fig.Q] for comparison, (b) I(V) characteristics of the 
UPd2Al3 PC at different temperatures, calculated according to Eq. © for d=200 nm 
and assuming p = 0. (c) Modification of the calculated I(V) at 1 K by adding the 
residual resistivity po to p(T). (d) Calculated /(V) curves at different temperatures 
for po = 10/^cm and d=200nm. 



the center of the PC depends on the applied voltage via [3J H] 

T 2 = TU + ^, (5) 

When the temperature T bu i k of the bulk sample vanishes, the contact temperature 
varies linearly with bias voltage like T = V/2yL. Using the standard Lorenz number 
L = L = 2.45 ■ 10~ 8 V 2 K~ 2 , a lmV bias voltage will raise the temperature of the 
contact by 3.2 K. 

In the thermal regime the /(V) characteristic of the contact depends on the 
temperature-dependent electrical resistivity p(T) according to |U E] 

I(V) = Vd f 1 dX (6) 

where T is defined by Eq. (j3J). We used the smooth curve in Fig. E[a) to approximate 
the experimental p(T), but omitted the SC transition. 
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Figure 6. Comparison between the experimental I(V) characteristics (solid curves) 
from Fig. compared to the calculated ones (dashed curves) from Fig. Efd) at low and 
at high temperatures. The bottom calculated curve is multiplied by 0.9 along the / 
axis. 

The calculated /( V) curves in Fig. EJb) had maxima at around 2 — 3 mV, which 
results in a hysteresis for up- and downward sweeps when the junction is driven by a 
current source. These maxima are the larger the steeper the drop in p(T) around T^. 
They decrease and become broader with increasing residual resistivity, see Fig. Efc). 

With voltage biasing we would expect to recover the full I(V) chracteristics without 
hysteresis. However, this would require to install small resistors near the sample in 
parallel and in series with the break junction to cut off its bistability, see for example 
Ref. ^7]. This was not practical in our experiments because in each cool down we 
wanted to investigate many break junctions over a wide range of resistances. 

Figure |H1 shows that the theoretical /(V) describe well the experimental data, 
including the width of the hysteretic features, using d = 200 nm and po — 10 pQcm. 
These are the only two adjustable parameters. And they fit well those parameters 
that have been derived above independently from the measured R(T). This agreement 
strongly supports our interpretation that local (at the PC) thermal effects determine 
the behavior of our UPd2Al 3 break junctions. 

Surprisingly, the standard Lorenz number Lq yields the best fit to the experimental 
/(V) curves while in bulk UPd 2 Al 3 L(T) rises from 0.6L below IK up to « 15 L at 
24 K because of the dominant heat transport by phonons ^Hj- This implies that at the 
PC the phonon channel is closed, and heat is carried away by electrons only. 

With increasing residual resistivity the hysteresis of the experimental /( V) curves in 
Fig. El transforms into an inflection point, corresponding to the maxima in dV/dl (inset 
of Fig. This trend agrees with the theoretical curves in Fig. EJc). However, it seems 
that with increasing po the experimentally observed I(V) hysteresis disappears more 
quickly than expected from theory. One could speculate that the larger p the more 
strongly degraded is the contact structure, simultaneously broadening the steep rise of 
p(T) around the AFM transition, see inset in Fig. 1 in Ref. [IB]. Such a broadening 
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Figure 7. I(V) characteristics of different UPC12AI3 break junctions at T ~0.1K. 
The current was swept upward (solid lines) and downward (dotted lines) . The residual 
resistivity po and the contact size d, derived from R(T), are given for each contact. 
Inset shows dV/dl for the bottom I(V) curve. 

would be similar to that of the SC transition of the PC in Fig. El 

The slightly reduced size of the hysteresis loops in an applied magnetic field, shown 
in the lower inset of Fig. EJ goes in the same direction. This could be attributed to a 
small positive magnetoresistivity of UPd 2 Al 3 [Tj5] . 

Note that the AFM transition itself is difficult to resolve in the I(V) characteristics. 
This transition shows up as a small step in the derivative of p(T). Since I(V) is described 
by an integral containing p(T) over a certain range of temperatures defined by the bias 
voltage, one would at least have to check the second derivative d 2 I/dV 2 . Nevertheless, 
the huge anomalies in I(V) reflect the AFM transition because the magnetic ordering 
improves dramatically the coherence of the electron scattering processes, leading to the 
steeply decreasing resistivity. 

In the /(V) characteristics SC appears as an 'excess' current, see the upper inset 
of Fig. El To calculate I(V) of the superconducting anomaly we assumed that p(T) 
varies like the contact resistance R(T) (see Fig-EJ) , normalized to the normal-state p . 
This lead to a single peak at around 0.15 mV, while the experimental /(V) in the upper 
inset of Fig. |3] rises almost continuously. Thus the thermal model jH Ej, developed 
for normal-state contacts, fails to describe even qualitatively the resistance of the SC 
contacts. This failure could have two reasons. First, the broad SC transition indicates 
that the contact does not have a single T c , but a whole distribution ranging from T c m 
at the center of the contact where it is reduced due to stress and disorder and T c = 1.8 K 
far away in the indisturbed bulk material. 

A multiply connected contact, where each connection has its own T c creating a single 
spike in the spectrum, could be excluded because Sharvin's resistance was recovered. 
Consequently also the normal-state residual resistivity may vary locally, both along 
and vertical to the contact direction, while Eqs. © and (jHJ) have been derived for 
homogeneous samples only. This would greatly affect the I(V) characteristic at low 
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bias voltages since the pattern of current flow could change abruptly, for example when 
the critical supercurrent is exceeded in part of the contact region. It will not change 
/(V) at large bias voltages because then the large T-dependent part of the electrical 
resistivity takes over. Second, at low temperatures, when both the elastic and the 
ineleastic electron mean free paths are largest, the UPd2Al3 contacts could be in the 
diffusive instead of the thermal regime. With increasing temperature or bias voltage the 
mean free paths get smaller, and the contact is forced into the thermal regime again. 

Break junctions with URu2Si2, another heavy-fermion SC, had peaks in the 
differential resistance dV/dI(V) at voltages described by Eq. (jSJ) [TU]. This indicated 
the destruction of SC in the constriction due to local heating. We have observed the 
same behavior also for our UPd2Al 3 contacts in the SC state. 

The UPd2Al3 junctions presented here are non- linear devices. Their iV-shaped /(V) 
characteristics have a negative differential resistance. Those devices could be applied - 
in principle - like an Esaki tunnel diode or a Gunn diode as amplifiers, generators, or 
switching units [213 E]. Of practical interest is therefore the possible minimum response 
time. We estimate it by the thermal relaxation time r ~ (cD/X)d 2 of the contact [3]. 
Here c is the thermal heat capacity per volume, D is material density and A the thermal 
conductivity. With the molar heat capacity of 3.5 J/K _1 mol _1 D ~ 10 gem -3 , and 
A ~ 4WK~ 1 m~ 1 J3| at 10 K the relaxation time becomes r ~ 100 ps for a d = 100 nm 
wide contact. This is three orders of magnitude larger than for a standard tunnel diode, 
but it could be reduced by using smaller contacts as long as they remain in the thermal 
regime. One (dis) advantage, however, is the low 4mV working point (at the maximum 
negative slope of I(V)), an order of magnitude below that of typical Esaki tunnel diodes. 

5. Conclusion 

Sub-/im scale metallic break-junctions of heavy-fermion UPd2Al 3 showed hysteretic 
I{V) characteristics at low temperatures. These highly nonlinear I(V) curves can be 
reproduced theoretically by assuming the constrictions are in the thermal regime. Such 
anomalous I(V) curves are due to the unusual p(T) dependence of UPd2Al 3 at the AFM 
transition. Since those point contacts with iV-shaped I(V) characteristics are nonlinear 
elements with a negative differential resistance, they could serve as the analogue of Esaki 
tunnel diodes or Gunn diodes as amplifiers, generators, and switching units. From this 
point of view UPd2Al 3 is not such a unique material - each metal with similar p(T) 
should produce similar /(V) characteristics. This can be expected for many materials 
which order magnetically at low temperatures, since their resistivity typically increases 
steeply when magnetic order is destroyed by thermal fluctuations. 
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